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Abstract Bondi 's approach to the construction of a coordinate system is used with a different choice of 
gauge, in accordance with which the radial coordinate r is an affine parameter, to cast the metric tensor 
into a form suitable for use with the Newman-Penrose null tetrad formalism. The choice of tetrad has 
the result that the equations and all the functions that appear in them are real-valued. A group 
classification of the Sachs equations in this gauge leads to a unique expression for the first of the five 
independent elements W of the Weyl spinor, and to the corresponding exact solutions for two of the 
metric functions on an initial null hypersurface. A proof is presented that the result for % constitutes 
the appropriate characteristic initial value function for all physically realistic axisymmetric, non-rotating 
vacuum space-times. Integration of the field equations on the axis of symmetry when an equatorial 
symmetry plane is also present, and on the equatorial plane itself when time rates of change can be 
neglected, shows that these data produce results consistent with Newton's laws and with the 
Schwarzschild solution in the appropriate limits. The solution on the axis of symmetry indicates that the 
Weyl curvature increases without limit between black holes as their separation decreases. 

1. Introduction 

It is well known that Einstein's vacuum field equations can be solved in principle for the metric tensor of 
a non-rotating axisymmetric vacuum space-time provided that certain initial data are specified. The 
problem of finding physically interesting solutions therefore reduces to that of finding the form of the 
appropriate initial data. It is the purpose of this paper to demonstrate that by combining Bondi's 
approach to the construction of a coordinate system [1],[2] with the Newman-Penrose (hereafter referred 
to as NP) spinor formulation of the field equations [3] it is possible to cast the equations into a particularly 
transparent form. With the use of a Lie group analysis an initial value function can be obtained from two 
of these equations. It will be shown (in appendix A) that this function provides the appropriate initial 
data for all physically realistic vacuum space-times that possess an axis of symmetry and are non-rotating 
(i.e. that have reflection symmetry with respect to rotation about the axis). One form of the initial value 
function yields the Weyl tensor that is implied by Newton's laws as an approximation on the axis of 
symmetry. It yields a pair of Schwarzschild solutions, and a second form yields a single Schwarzschild 
solution, as an approximation in the appropriate case. 



2. Coordinate system and tetrad structure. 

Bondi's coordinates are based on a set of null cones generated from the timelike curve followed by a point 
O that remains on the axis of symmetry. The null directions at O are labelled with spherical coordinates 
9 and <|). Radial distances are measured by the parameter r that characterises individual generators of the 
null cone, and the coordinate set is completed with a retarded time v, which is constant on each null 
hypersurface. Assign numerical references 1,2,3,4 to v, r, 8, and § in that order. "Geometrised" units of 
measure, in which c = G = 1 , and a space-time signature of -2, will be employed throughout. 

Null geodesies are characterised by v, 8, <|) = const. It follows from the fact that only r varies along an 
outgoing null geodesic that g 2 2 must be zero. Symmetry under the reflection § -> -§ requires that 

El4 = §24 = §34 = 0. 

The absolute derivative of the tangent to a curve, along the same curve, reduces to the following 
expression if the curve parameter r is used as one of the coordinates: 

i»v/ v = r v 22 (2.1) 
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It follows that if the curve is a geodesic r 2 2 = r 22 = • Bondi et. al. [2] demonstrate that this 
implies that the ratio g23/gi2 must be independent of r. This allows g 23 to be eliminated by the 



transformation: 

( 

(2.2) 



du = X 

where A is a suitable function of v and 6. 



dv + ^1 de 



§12 



Bondi et. al. use the remaining coordinate freedom to make r an area distance. If the coordinate freedom 

is used, instead, to transform r in such a way that g i2 = 1, the consequence is that r 22 = and r 

becomes an affine parameter. The metric tensor now has the canonical form used by Newman, Penrose 
and their co-workers. These coordinates look like Robinson- Trautman [4] coordinates but there is a 
significant difference in the fact that the current coordinates are related through (2.2) to a set of null cones 
that have their vertices on the axis of symmetry, whereas Robinson-Trautman coordinates do not require 
the existence of null cones. As pointed out by Penrose [5], the presence of the vertex on a null cone 
provides a common reference point to which field variables may be referred by parallel transport along the 
generators of the null cone. 

If it is assumed that the manifold has a spinor structure [6], a null tetrad may now be introduced at every 

point of the space-time manifold, using as reference the tangent vector to the null geodesies that form 
the coordinate rays. A spinor structure requires that the manifold should have a unique orientation at 
every point. To achieve this, the neighbourhood of the origin must be simply connected and may not be 
extended beyond any points at which the coordinate rays form cusps (i.e. where the null geodesies cross 
one another). This means that in general the space-time will have disjoint coordinate patches. 

The nature of the coordinates and the form of the metric tensor established above means that: 

^ = 8£, ^ = 5^ (2.3) 

where 8 is the Kronecker delta. 

The tetrad is completed with a second real null vector rt l normalised so that n^ = 1 and a complex 

null vector nv 11 together with its complex conjugate m^ (for the details see NP [3] ). The tetrad vectors 
are related to the elements of the metric tensor by the relation 

A similar relation holds for the components of the covariant form of the metric tensor. 
To satisfy the orthogonality relations the tetrad vectors must have the following form: 

n» = 8J 1 + U82 1 + (i = 3,4) (2.5) 

and m" = + $8? (2.6) 

where, using the notation of NP and their co-workers, U, X 1 , co , and E, 1 are unknown, complex- valued 

functions of u,r, and 6. Using (2.4) together with (2.5) and (2.6) gives the elements of the metric tensor 
in terms of these functions. 

The vector can be chosen so that at the origin its spacelike projection lies in a plane containing the axis 
of symmetry. Symmetry then requires that it must remain in this plane as the tetrad propagates outwards. 

One can therefore put X 4 = so that there is no ambiguity in putting X 3 = X . The complex null 
vector rrf is defined relative to a pair of real vectors which span the spacelike hyperplane tangent to the 
surfaces u = const. If the imaginary axis is chosen to lie along the direction of the coordinate x 4 then 
^ 3 becomes a real-valued function. It is easy to show that all the metric functions are then real-valued 

except for ^ 4 which is pure imaginary . 

For later convenience put = 1/q and E, 4 = i/Q where both q and Q are real- valued functions of the 
coordinates. The metric tensor resulting from this analysis is: 
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ds z = (-2U - \ X 2 q 2 + 2Xcoq) du z + 2 du dr - (2coq - q z X) du d9 - | q z d9 z - ± Q z d<|> 



1 ^2 ja2 1 r.2 2 



1 2(U-co 2 ) X-^ 



X 



2co 








Q 2 ^ 



(2.7) 



As is the case with Bondi's metric tensor this metric contains essentially four unknown functions. The 
apparent redundancy following from the fact that five unknown functions survive from the specification of 
the tetrad vectors will be seen later to disappear when the condition for parallel propagation of the tetrad 
is imposed. 



It is easy to see that the functions co and X measure the focusing of the coordinate rays. For an 
asymptotically flat space-time they must tend to zero at infinity. When it can be assumed that the metric 
tensor will approach the Schwarzschild form as r -> go , then the affine parameter r must at the same 
time tend to the Schwarzschild radial coordinate (which is an affine parameter for radial null geodesies in 
the Schwarzschild space-time). Also, in this case one must have 

U -»-- + — , q^V2r, Q^V2rsin6 as r -> go (2.8) 
2 r 

where m tends to the Schwarzschild mass of the system. 



3. The Newman-Pen rose equations 

Expressed in terms of the metric tensor (2.7), all the spin coefficients, all the elements of the Weyl spinor, 
and all the NP equations are real-valued. This result greatly simplifies the equations without any loss of 
transparency. The eight "metric equations" that relate the spin coefficients to the elements of the metric 
tensor are all satisfied identically. To avoid repeating the numerous NP equations themselves they will be 
referred to by the equation numbers employed in [3]. 

The condition for parallel propagation of the tetrad is that the spin coefficient n must vanish. This 
provides the relation: 

qSX + »Sq + So^ Q ^ 
2 9r q dr dr 

The NP equations (6.1 la) to (6.1 lr) are reduced to 12 equations when (3.1) is taken into account. It can 
be demonstrated that any solution of these equations, in the current coordinate system, will be regular on 
the axis of symmetry. This demonstrates consistency with the assumption of regularity inherent in their 
derivation and that of the metric tensor (2.7). The analysis also produces a number of results that are 
required in section 5. 

The integrability of the equations on the axis of symmetry follows trivially from Frobenius's theorem [7]. 
(Note that, while there is a coordinate singularity on the axis, points on the axis can be included in the 
u,r,6 hyperplane by analytic continuation.). The behaviour of the metric tensor and the tetrad vectors on 
opposite sides of the axis of symmetry shows that the functions co, X, and Q must be odd under the 
reflection 6 -> -6 . The definitions of the elements of the Weyl spinor show that the same must be true 
of the functions *Pi and ¥3 . If these functions are analytic then the functions themselves and all their 
even partial derivatives with respect to 6 must vanish on the axis. The opposite is true of the functions q, 
U, *P , ^2 , and ¥4 . These functions must be even so that all their odd partial derivatives with respect to 
6 must vanish on the axis. 



The spin coefficients p and a, also called the "optical scalars", provide additional information. In the 
current coordinates the definitions of these functions lead to the following expressions: 

P + o = -i§ " " (3.2) 

q or 
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Sachs [8] first showed that o measures the shear of the null congruence while p measures the 
expansion of the congruence. On the axis of symmetry o must be zero. 

Equations (6.1 la) and (6.1 lb) in NP were discovered by Sachs [8] . They read: 

f = p 2 + a 2 (3.4) 
or 

^ = 2pa + ^ (3.5) 
Sr u 

It follows immediately that on the axis of symmetry 

Y = and P = ^ 1 — (3.6) 
P -r 

where p° is a function that is independent of r. The first of these two results also follows directly from 
the Petrov classification of the space-time on the axis of symmetry*. The possibility that p = 
everywhere can be excluded immediately on the basis of elementary flatness, since a null cone in flat 
space is characterised by p ^ . 



Subtracting (6.1 Id) from (6.1 le) in NP, multiplying by Q and then allowing the functions to tend to their 
values on the axis of symmetry gives a first order ordinary differential equation for 9Q/96 (regarded as a 
function of u and r defined on the axis of symmetry). Without loss of generality the analysis can be 
confined to the first quadrant so that the value of 6 on the axis is single-valued and constant, with the 
value 8 = 0. Solving the equation then gives the result that 9Q/96 = f(u) q where f(u) is a function of 
integration. Multiplying (6.11m) by Q and letting the functions tend to their values on the axis gives 
df/du - 0. It is possible therefore to adopt f(u) = 1 as an initial value, so that for all u on the axis of 
symmetry: 




(3.7) 



This is the condition for regularity on the axis of symmetry. Its validity makes it possible to use 
l'Hopital's rule to evaluate the indeterminate forms 0/0 that appear in the NP equations when Q is in the 
denominator and the equations are evaluated on the axis in the current coordinates. 



4. Null hypersurface initial data 

The NP equations can be divided into 'hypersurface equations' that involve no derivatives with respect to 
u, and the remaining equations that do involve u derivatives. The equations can be integrated by starting 
from an initial function, specified on a single null hypersurface. The values of all the metric functions on 
the initial hypersurface can be determined from the hypersurface equations. The first derivatives with 
respect to u can then be determined from the remaining equations. Higher u derivatives can be 
determined by differentiating all the equations with respect to u and repeating the cycle. This allows a 
formal expansion in powers of u to be generated for each of the metric functions. The conditions under 
which such expansions converge have been investigated by Friedrich [14]. 

In the current formulation of the field equations the natural candidate for an initial function is % , the first 
of the five independent elements of the Weyl spinor. This follows from an analysis carried out by 
Penrose [5] who demonstrated that the values of *P on the points of an initial null hypersurface in an 
analytic empty space-time, together with Einstein's field equations, determine the entire field. Penrose 
also demonstrated that these initial data can be chosen freely without any constraints (when there are no 
constraints on the geometry of the space-time), except perhaps that they should be analytic. 
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The hypersurface equations that provide the first step in the integration of the field equations in the current 
form are the Sachs equations (3.4) and (3.5). In the coordinates used here they take a particularly simple 
form: 

f^ + Y q = and ^-^ Q = (4.1) 

Since *F is a function of r and 6 that is not determined by the NP equations, it can be specified 
independently of these equations on the initial null hypersurface. Also, the two Sachs equations involve 
differentiation with respect to r only. They can therefore be treated as a closed system of ordinary 
differential equations that determine the r-dependence of q and Q once the character of the space-time has 
been prescribed by the choice of % on the initial hypersurface. Much attention has been directed at 
determining the nature of asymptotic expansions that constitute physically realistic choices for the form of 
the initial data, mainly with a view to exploring the nature of gravitational radiation at infinity. NP and 
their co-workers chose an expansion in terms of powers of 1/r with the condition that % = O (r 5 ) for a 
general (not necessarily axially symmetric) radiation metric. More recently it has been argued [9] that 
polyhomogeneous space-times, admitting expansions in terms of r" J log'r , may be more appropriate. 



As an alternative to prescribing a series expansion of *P it seems natural to carry out a group 
classification [10] of the Sachs equations (4.1) to determine whether there are any functional forms for % 
that would extend the principal Lie algebra possessed by the pair of equations. (The principal Lie algebra 
of a set of differential equations is the Lie algebra possessed by the equations regardless of the form of an 
undetermined function.) On the initial null hypersurface the value of % is independent of q, Q, and their 
derivatives [5]. It is therefore possible to put % = % (r,0) in the space of differential functions [10] for 
these equations. It follows that the principal Lie algebra of (4.1) corresponds to the pair of scale 
transformations of q and Q. This is consistent with the fact that Einstein's field equations are invariant 
under scale transformations of the elements of the metric tensor [11]. (However, note that this holds for 
the full set of field equations only when the same re-scaling factor is applied to all the elements of the 
metric tensor.) It turns out that if % is regarded as invariant, additional symmetries appear if, and only 
if, *P takes the following form [12]. The value u = designates the single null hypersurface on which 
the result applies. 

y ol„=o= rr 2 [° - \ 2 (42) 
IV + V + bo) 

where a and the bj are functions that are independent of r. 



The equations then possess a three-dimensional Lie algebra spanned by the operators 

2 . , ,3 d ^ d 
X!=(b 2 r +b 1 r + b )— + rq— + rQ — 
or dq dQ 

3 a (4J) 

where X 2 and X 3 constitute the principal Lie algebra. 

The analysis in appendix A demonstrates that the hypersurface equations for any solution of the field 
equations representing a physically realistic non-rotating axisymmetric space-time must be invariant under 
a Lie point transformation that transforms the coordinate r as the operator X] does, and that leaves the 
function % invariant. The only functional form of % that allows (4.1) to admit such a symmetry is (4.2). 
It follows that % must take the form (4.2) on an initial null hypersurface for any solution representing a 
physically realistic non-rotating axisymmetric vacuum space-time. The analysis in appendix A also 
indicates that the bj must be constants. The function ao arises as a function of integration and is, in 
general, a function of 6. 

Penrose's result [5] concerning initial null data is extended by a theorem due to Friedrich [14] that proves 
the existence of analytic solutions of Einstein's field equations for analytic initial data on past null infinity 
and on an incoming null hypersurface. (Note that the vertex corresponding to any particular null cone is a 
point that is, strictly speaking, not part of the null hypersurface associated with that cone. Regularity 
conditions therefore have to be imposed at the vertex to reflect, by analytic continuation, the regularity 
that is assumed to exist on the null hypersurface associated with the past null cone.) It follows that there 
are exact analytic solutions of Einstein's field equations corresponding to the initial null data provided by 
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(4.2). These solutions can be found, in principle, by the procedure outlined above, but in practice the 
complete integration of the equations remains a formidable computational problem. 



The solutions of (4.1) with % in the form (4.2) take five different forms depending on the relative 
magnitudes of ao and the bj . Assuming b 2 * then it is possible to set b 2 = 1. If 



4b + 4 a < the solution can be put in the form: 



q°(r 2 + bjr + b Q )z cos 



4(b Q +a Q )-b; 
4b n -b? 



arctan 



4b 







2r + b 



(4.4) 



where q and <|) q are functions of integration that are independent of r. 



The corresponding solution for Q is similar except that the sign of a is reversed. For the case 

bj - 4b Q b 2 — 4 1 a Q | > the polynomial in the denominator of (4.2) has two unequal roots at which *F has 

singularities. Putting b 2 = -1 (so that the expression under the square root sign remains positive in the 
application to be discussed below) , the solution in this case can be put in the form: 



= q°(b +b 1 r-r 2 ) 2 sinh 



4 bf + 4b 



ln^ 



OJ 



(b 2 + 4b ) 



+ 4b n 2 +b, -2r 



21 



4(b +b ir -r 2 ) 



(4.5) 



V I J J 

The corresponding solution for Q is similar except that the sign of ao is reversed. Other forms of the 
solution are given in appendix B. 

When the coordinate r can extend to infinity without encountering singularities, as is the case when the 
solution takes the form (4.4), then q and Q tend to their flat-space values as r tends to infinity while ao and 
the bi remain finite (assuming appropriate values for the functions of integration). These are also the 
values q and Q take for the Schwarzschild solution. It is easy to show that the Schwarzschild metric is 
indeed a trivial solution of the equations corresponding to the case a = 0. 



5. The exact solution on the axis in the initial null hypersurface when bi = and there is 
reflection symmetry about the origin with respect to r. 

The transformations corresponding to (4.3) are one-parameter diffeomorphisms. Accordingly, they 
describe coordinate transformations within a single coordinate patch. In the case of a pair of black holes 
moving along the axis of symmetry, focusing of the coordinate rays would give rise to a coordinate patch 
that encloses an interval of the axis itself between the pair but does not extend to infinity. The solution on 
that part of the axis enclosed between the two black holes would therefore only describe the gravitational 
field measured by observers who are themselves on that part of the axis. 

When the origin lies between a pair of black holes, there will be singularities on the axis on both sides of 
the origin. It is natural therefore to use (4.5) for this case. It follows from (4.2) and (3.6) that ao is zero 

on the axis of symmetry. When ao tends to zero, (4.5) can be brought into the form q = V2r by 

assuming that bi = 0, that <|) q tends to zero as 6 -> , and that q° tends to the value V2 . It follows 
from (3.6) and (3.2) that this value of q on the axis is appropriate to an observer who coincides with the 
point O from which the null cones that define the coordinate system are generated 

With q independent of u on the axis of symmetry, the NP equations can be solved for 9X/96 and 
dm/dQ on the axis using NP equations (6.1 If), (6.1 lg) + (6.11h) and (3.1) above. The resulting 
expressions can be used in (6.1 11) (after making use of l'Hopital's rule) combined with (6.1 If) to obtain a 
second-order linear differential equation for U. The solution on the axis of symmetry is: 
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Til - 1 f 5 Q 5 q dr r f 1 



1 [ S 3 Q 5 2 q 



dr + ^f + ^ 2 - + k 3 r (5.1) 
x 1 r 



where the kj are constants of integration. The integrands can be evaluated using (4.5), taking account of 
the symmetry properties, discussed in section 3, of Q and q and their derivatives with respect to 6 and 
of (3.7). 

Assume b is non-negative so that one can put b = d 2 . The constants of integration ^ and k 2 can be 
evaluated using l'Hopital's rule and the condition that U must have a finite limit as r -» . The function 
A in (2.2) can be normalised so that it is unity on the axis. It then follows from the construction of the 
coordinate system that U = - 1/2 atr = 0. Using this condition it follows that: 

56 56 d 56 

If one assumes reflection symmetry with respect to r about the origin one can require that SU / Sr -» as 
r -» . This condition allows k 3 to be determined. 

From (6. 1 1 f) in NP, and the fact that co must be zero on the axis of symmetry, it follows that: 

Using the above results in (5.1) and then substituting in (5.3) gives an expression for ^2 that is well- 
behaved as r -> and gives the limit: 

1 52a ° 1 (5.4) 



lim *¥ 2 

r^O 



=0 6 d ae 2 d 3 

This expression can be compared with the Newtonian equivalent for a pair of identical massive bodies, 
each of mass m, located at positions +s and -s relative to an origin at their centre of mass. This can be 
calculated from the geodesic deviation equation applied to a pair of test particles separated by an 
infinitesimal spacelike displacement along the axis of symmetry and initially at rest at a distance r from 
the centre of mass. Taking into account that in an empty space-time the Weyl tensor equals the Riemann 

tensor, and that from the definition of in the tetrad employed here *F 2 = ~^C 1212 on the axis of 

symmetry, one finds, since the four-velocities of the test particles can be considered to be orthogonal to 
the axis of symmetry at r = 0: 

lim *F | (Newton) = -mi — l — + — l — (5.5) 
r^O 2 ' e - U [(s + r) 3 (s-r) 3 J 

If d is identified with s and one puts: 

5 2 a 

^- = -12md (5.6) 

56 

then the limit (5.4) is identical to that of the Newtonian expression (5.5), at r = 0, on the null hypersurface 
under consideration. This is a satisfying result, considering that ao must be zero when m = and when 
d = . Note also from (5.2) that, with this result, the 6 derivatives of q° and Q° tend to their flat space 
values as m/d -> . 

Substituting for the kj in (5.1) and using the identification (5.6) yields the metric function U on the axis of 
symmetry in a null hypersurface for the maximally symmetric case dealt with here, i.e. that of a pair of 
identical non-rotating black holes located on an axis of symmetry with the origin at their centre of mass: 

1 2m/ ? „ ,2V ( d + r^| 4md, a, w, n 8mdln(d) 8m 1 .. _. 

U u ,e=o=— 2 "3d 2 In — l n {(d + r)(d-r)} + (5.7) 

rd 2 x 1 x 1 d 2 

The value of U corresponding to the Newtonian approximation (5.5) for ^2 can be obtained by using 

(5.3). Choosing the constants of integration appropriately, this yields : 

1 y m m 2m 1 

U L_n (Newton) = + (5.8) 

l0 -° d + r d-r d 2 
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This clearly approximates to a pair of Schwarzschild solutions, each of mass m, centred at r = ± d , for d 

large compared with m and d - r small compared with d + r. Taylor expansions about the value r = 
show that (5.7) and (5.8) are identical for small r/d, differing only from the fourth order terms onwards. 
This means that the exact relativistic metric function (5.7) would be indistinguishable from (5.8) in a 
sufficiently small neighbourhood of the origin. 

Note that (5.7) tends to the value 8m(l-ln2)/d-^ as r -» . The fact that U does not become zero 
for r < d until d » 4.9 m can be interpreted as follows. For a coordinate ray at an infinitesimal but 

non-zero angle 6 the metric function q is given by (4.5). With b! zero and b = d 2 , q will have two zeros 
at positions r = ± d . This means that the null geodesic congruence emanating from the origin has cusps 

at these points. The fact that U, as given by (5.7), is greater than zero until relatively small values of d 
can be interpreted to mean that for larger separations of the pair of black holes a cusp is encountered along 
the axis before a point of infinite red shift is reached. 

Substituting from (5.6) in (5.4): 




(5.9) 



This expression is identical to the Newtonian expression for a pair of point masses each of mass m, 
separated by a distance 2d. However, note that in the relativistic case 2d is measured in terms of an affine 
parameter, and is the distance that separates the pair of points on a null geodesic along the axis of 
symmetry that are conjugate to the origin. An interesting conclusion follows from this result, since the 
null hypersurface on which the result holds is arbitrary: the "tidal force" component of the Weyl tensor on 
the axis of symmetry between a pair of non-rotating black holes in head-on collision increases without 
limit as the points on the axis at which there are cusps in the null congruence emanating from the origin 
approach one another. This result is probably valid more generally than for the maximally symmetric 
case considered here 



6. The solution on the equatorial plane when there is reflection symmetry about the origin and 
time rates of change can be neglected. 

The symmetry relations on the equatorial plane, which is a plane of reflection symmetry in the case being 
considered, are: 

Sq SQ 5U 

X = co = — - = = 

56 







(6.1) 



56 56 

A second order differential equation for U in terms of q and Q can be obtained from NP equations (6.1 lg) 
- (6.1 lh) by substituting for W 2 from (6.1 If) and by using the symmetry relations in (6.1). Assuming 
time rates of change can be neglected the solution is exactly Schwarzschild when q and Q take the forms 
(2.8). Using the non-singular forms (4.4) for q and Q the solution is: 

l ( \\ 



U 



-1 — 

1 Q 



c 2 q 



tan < 



Q(4a + 4b 







1 + 



4 a, 







4b 







arctan 



h2 



(4b -bf) 
bj + 2r 



(6.2) 



where C! and c 2 are constants of integration and the symbol 1 indicates that the equation applies on the 
equatorial plane in the initial null hypersurface. This expression can be expanded as a Taylor series in 1/r 

= _ l 
2 

2mb, 



around the point 1/r = 0. Choosing Ci 



and c 2 = -4m the result is: 



U 



1 2m 



*0 



8ma + 3b ] a 



3mb 



2 



(6.3) 



2 r r 2 r J 

To first order this is a single Schwarzschild solution with mass 2m when the functions that are 
independent of r are finite. 



7. Conclusion 

Although the complete solution has been computed only in the initial null hypersurface for the simplest 
cases of the general axisymmetric problem, the results clearly show that the characteristic initial value 
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function (4.2) produces expressions that are compatible in the appropriate limits with Newton's laws and 
with the Schwarzschild solution. To progress beyond the simplest cases, substantial computing power is 
required and it is necessary to evaluate difficult quadratures. 

The prospect of detectors capable of measuring gravitational radiation from astrophysical events has led to 
much work on numerical solutions of Einstein's equations [15]. This work is ultimately directed at 
understanding the gravitational radiation expected to be produced by a pair of black holes spiralling in 
towards one another. Work on vacuum solutions completed to date has included the case of a pair of 
non-rotating black holes in head-on collision, starting from initial data that represent "two nonrotating 
black holes poised near one another in a momentarily stationary configuration" [16]. These initial data 
are based on an exact solution on a hypersurface of time symmetry, obtained by Misner [17]. The 
solution to the initial value problem reported here relaxes the requirement of time symmetry. 

Some workers have used approximate initial data [18] for numerical solutions. However, the uncertainty 
inherent in inexact methods of solving Einstein's equations, and the importance of obtaining exact 
solutions wherever possible, has been stressed by Ellis [19]. An exact analytic radiation solution can be 
obtained, in principle, from an extension of the work reported here, as outlined in section 4. If this can be 
achieved in practice, even for a restricted spatial geometry, such as that of section 6, the result will serve 
as a validation test for numerical codes. 

Numerical methods are inherently not suited to the investigation of fields in the vicinity of singularities 
and the curvature singularity revealed by expression (5.9) has not previously been reported. Recent 
evidence concerning the presence of black holes in galactic nuclei [20] and evidence linking galactic 
collisions to quasars [21] suggest that the exact solution for the strong and rapidly varying fields in the 
vicinity of this singularity is likely to lead to valuable insights concerning physically observable processes. 
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Appendix A: Symmetry transformations of the hypersurface equations. 

The particular forms of equations (4.1), together with those corresponding to the remaining NP vacuum 
equations expressed in the coordinates employed in this paper, are determined entirely by the assumptions 
used in section 2 to construct the coordinate system. These assumptions are not affected by the location 
from which the null cones, which define the coordinates, are observed. A surface v = const can be 
marked by a physical light pulse. Provided the metric of the space-time is known, observers positioned at 
different locations can follow the surfaces v = const as they expand through space by monitoring the 
reflection of the corresponding light pulses from test objects. Their measurements will be consistent with 
the NP equations. Under the following conditions the transformations that map the measurements of the 
observers onto one another will constitute Lie point symmetries of the NP 'hypersurface' equations. 

Let the observers all be members of a set S who are located on the axis of symmetry in a regular 
neighbourhood of O and let them use clocks that have been synchronised at a single event when they all 
coincided. The axis of symmetry defines an invariant direction common to all the members of the set S. 
Since g 2 3 is zero on the axis of symmetry for every member of the set S in the original coordinates, before 
applying the transformation (2.2), the same physical light pulse v = const, can be used to identify and label 
the corresponding hypersurface u = const in the frame of reference of each observer. If, furthermore, 
each observer uses a boost to match the intervals between light pulses, emitted at intervals dv and received 
directly from O, to his own clock, then each observer can set the function X in (2,2) to unity on the axis, so 
that du = dv for every observer in the set S. Within the single null hypersurface under consideration the 
different observers are then distinguished only by a single parameter, being their relative separation along 
the axis. 

An observer will, in general, assign the label u = const to a pulse of light emitted at a different time for 
each angle, such that the intersection of the null surface u = const with a surface of simultaneity, or 
spacelike slice, in the coordinates of that observer (the 'wave front'), constitutes a sphere on which the 
coordinate r is the same on the null geodesies corresponding to every angle. This follows from the fact 
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that when g 23 is zero, the direction (9/96)^ is orthogonal to the direction (9/9r)^ . In the coordinates 

of each observer, the measurements of r, 6, and § on these surfaces u = const will be related to the values 
of q, Q, and % on the same surfaces by (4.1). 

Transformations from the reference frame of one observer in the set S to that of another can be made as 
follows. The value of r can be deduced from measurements carried out on the wave front v = const in a 
small neighbourhood of the axis of symmetry. Since g 23 is zero on the axis of symmetry for every 
observer in the set S it follows that for a spacelike displacement ds within the hypersurface v = const: 

lim ds 2 = --q 2 d6 2 (Al) 
e^o 2 

Since o = on the axis of symmetry, it follows from (3.6) and (3.2) (setting p° to zero) that 

limq=q°r (A2) 

where q° is a function that is independent of r. If one confines the argument to the first quadrant, as is 
possible without loss of generality, then 6 is single-valued and constant on the axis of symmetry. On the 
axis, q° can therefore only depend on u. It follows that: 

,. ds .... 
r = — lim — (A3) 

q ° e^o d6 

Points in the hypersurface can be located by specifying s and 6 as independent variables, instead of r and 
6, with s = 6 = on the axis. (There is a coordinate singularity on the axis when using s and 6, analogous 
to the coordinate singularity at the origin of a spherical coordinate system. However, values of r 
corresponding to events on the axis in the wave front can be assigned by analytic continuation.) Under a 
change of basis, in general, s' = s'(s,6) and 6' = 6'(s,6) . It follows that under a change of basis within 
the set S: 

9s' 9s' 
Jl ds' Jl — ds + — d6 

r' = — lim = — lim — (A4) 

q° d6' q° e^O 3e_ ds+ 96; de 

9s 96 

If the limits of the partial derivatives in the above expression are independent of r, depending only on the 
relative positions of the observers, then they are constants and the expression can be re-arranged to read: 

ctr + 6 

r' = (A5) 

yr + 8 

where a, (3, y, and 8 are constants. Without loss of generality the condition ct8 - y (3 = 1 can be imposed. 
This is then the three-dimensional Lie group SL(2,R). 

This group of transformations leaves the NP equations unchanged. However, the function % is not 
determined by the NP equations [5]. In general *P will be changed by the transformation (A5) so that the 
transformation, when applied to a solution of the NP equations in the current form, will yield a different 
solution of the same equations. This is a particular case of a general result, discovered by Geroch [13], 
that the action of the group SL(2,R) on any exact solution of Einstein's vacuum field equations with a 
Killing vector produces another exact solution that also possesses a Killing vector. Geroch's result 
confirms that the partial derivatives in (A4) must be independent of r. 

By construction, the map from the manifold defined by the set of points 2 = (r, 6 | u, § - const) , 
corresponding to any observer in the set S, to the manifold defined by the set of points 
2' = (r', 6' |u', <|> = const) corresponding to any other observer in the set S, is a diffeomorphism. (To 
give the sets 2 and 2' a manifold structure the singular point O must be excluded). Transformations 
from 2 to 2' therefore follow the normal tensor transformation laws. It follows that *P , which is a 
scalar function, is form- invariant under a transformation from 2 to 2' . 

This result means that both the NP equations and the form of *P on the initial null hypersurface must 
remain unchanged under the action of a transformation from the measurements of one observer in the set S 
to those of another. Since these transformations conform to the axioms for a group of Lie point 
symmetries, the corresponding transformations of the dependent variables q and Q must be part of a Lie 
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group. The result (4.3) shows that, to meet this requirement, the transformation of r must consist of a 
linear combination of the elements of the group SL(2,R), and also that the dependent variables q and Q 
have to transform, uniquely, as determined by the infinitesimal generator X] . This restricts the form of 
% to (4.2) and limits the set of solutions of the NP equations accordingly. 

Formal solutions of the NP equations may exist that correspond to functional forms of % different from 
(4.2), but they will not be physically realistic representations of the geometry of the space-time under 
consideration in this paper. The above analysis has therefore proved the following theorem: 

Theorem: Let M be a simply connected, analytic neighbourhood of a point O, on the axis of symmetry, 
that generates a timelike curve in an axially symmetric vacuum space-time in which Einstein's field 
equations are satisfied. Let the space-time possess reflection symmetry with respect to rotation about the 
axis of symmetry. Let M satisfy the axioms for a space with spinor structure. Consider a null 
hypersurface u = const constructed according to expression (2.2) from a sequence of null cones 
generated from O. Then, ifr is an affine parameter that measures radial distances from O along the null 
geodesies that pass through O, the r-dependence of the values of the component of the Weyl spinor on 
the null hypersurface, u=const, in M is given by the function (4.2) where the b t are constants and a is 
a function of 6. 



Appendix B 

The three cases of the solutions of the Sachs equations with % in the form (4.2) not mentioned in the text 
are as follows, putting b 2 = 1 : 

Ifb 1 2 -4b Q <0 but bj -4b Q +4|aJ >0,then 



Q = Q°(r 2 +b 1 r + bJ 2 cosh 



q = q° (r 2 + h^r + b Q ) 2 cos- 



b ,- 4 ( b o- a o) 



4b o" b i 



4 ( b o +a „)- b i 
4b -b 2 



arctan 



arctan 



( 4 W 
2r + h 



1 - 




( 4b o- b i 2 )' 


+ <b > 
q 


2r + bj 



if ao > 0. If a < the forms of q and Q are interchanged. 
If the quadratic has two equal roots, i.e. b\ = 4b then: 

V 2 1 

\ 



Q = Q | r + -± |cosh 



\ 



':;=q°|r + y |cos 



— r L + * 
r + J-b q 

V 2 1 



for a > 0. If ao < the forms of q and Q are interchanged and the sign of a is changed. 

2 I I 

If the quadratic has two unequal roots a and b and \> x - 4b Q - 4|a Q | < then 

l 



q = q [(r - a)(r - b)] 2 cos< 



(a-b)' 



Q = Q [(r-a)(r-b)]2 cosh< 



1 



2 , fr-a 
In, 

r-b 



A 



.(a-b)' 

if ao > 0. If a < the forms of q and Q are interchanged and the sign of ao is changed . 

In all cases cosines can be replaced with sines and vice versa. Also, since the equations are linear, a sum 
of the cosine and sine forms can be used as a solution in each case with the result that the hyperbolic cases 
can take the form of exponentials. 
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